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Abstract 



Harmonic generation by an atom in a laser field is described by the three-step 
mechanism as proceeding via above-threshold ionization (ATI) followed by the 
electron propagation in the laser-dressed continuum and the subsequent laser 
assisted recombination (LAR) . An amplitude of harmonic production is given 
by the coherent sum of contributions from different intermediate ATI channels 
labeled by the number m of absorbed laser photons. The range of m-values 
that gives substantial contribution is explored and found to be rather broad 
for high harmonic generation. The coherence effects are of crucial importance 
being responsible for the characteristic pattern of harmonic intensities with 
a plateau domain followed by a cutoff region. Due to multiphoton nature of 
the process, an efficient summation of m-contributions can be carried out in 
the framework of the saddle point method. The saddle points correspond to 
some complex- valued labels m = rric associated with the intermediate effective 
ATI channels in the three-step harmonic generation process. The advantage 
of this approach stems from the fact that summation over large number of 
conventional ATI m-channels is replaced by summation over small number of 
effective mc-channels. The equation governing rUc has a transparent physical 
meaning: the electron ejected from the atom on the first (ATI) stage should 
return to the core to make LAR possible. The effective channel labels m move 
along characteristic trajectories in the complex plane as the system param- 
eters vary. In the cutoff region of the harmonic spectrum a single effective 
channel contributes. For lower harmonics, in the plateau domain, two effec- 
tive ATI channels become essential. The interference of their contributions 
leads to oscillatory pattern in the harmonic generation rates. The calculated 
rates are in good agreement with the results obtained by other approaches. 

PACS numbers: 32.80.-t, 42.65.Ky, 32.80.Rm, 32.80.Wr 
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I. INTRODUCTION 



A. Three-step mechanism of HHG 

High efficiency of various processes in strong laser field could be understood basing on the 
observation that the field-induced quiver motion supplies an electron with high instantaneous 
energy. Rescattering of the energetic electron on atomic core generally is accompanied by 
the energy exchange between electron, core and electromagnetic field. In particular, the 
core could be excited or ionized [double ionization (DI) of an atom] and the high-energy 
photons could be emitted. The latter process is known as the harmonic generation (HG). If 
one considers the active electron initially bound to the core, then the electron at first should 
be (virtually) released to the quasi-free state, as a precondition that the subsequent events 
listed above become possible. As a whole, this constitutes three-step mechanism comprising 
above-threshold ionization (ATI), propagation in laser-dressed continuum and the final step 
which is electron-atom impact in case of rescattering, core excitation or ionization, or laser 
assisted recombination (LAR) in case of HG process. Strong interaction between the receding 
electron and the core is omitted in the standard Keldysh model of multiphoton ionization. 
Its importance was first pointed out by Kuchiev 0, who had predicted several phenomena 
where the electron-core interaction plays a crucial role. The related mechanism was named 
"atomic antenna" to stress the role of active electron in gaining energy from laser field. 

Specifically for HG the three-step model was promoted in the hybrid classical-quantum 
framework by Corkum ||^, see also the papers by Kulander et al [^,^. The subsequent 
theoretical developments were based on more sophisticated approaches and led to important 
advancements [|5HTT[], albeit the three-step nature of the HG process was somewhat veiled 
in these formulations. Most clearly the three-step mechanism is exposed in the framework 
of factorization method developed in Ref. |]T2|,ITE[ as an implementation of idea of atomic 



antenna of Ref. and applied to DI process. This technique allows one to present quantum 
amplitude as a direct sum over contributions of intermediate ATI channels. The application 
to HG process worked out by Kuchiev and Ostrovsky |[T^JT5[] quantitatively demonstrated 



validity and power of this approach. Since this development is the starting point of the 
present study, we briefly summarize it in Sec. below. 

The intermediate ATI channels form a discrete set being labeled by a number m of 
absorbed photons. Generally speaking, entire manifold of m-channels contributes to the 
rate of three-step process under consideration. One can anticipate that actually only some 
effective m-range should be essential, but this issue has not been investigated previously. 
Another important question arises from the fact that contributions to the rate of three- 
step process from different ATI channels are to be summed coherently. How essential are 
coherence effects in reality remained unclear. The acuteness of this problem is enhanced in 



view of recent publications by Becker and Faisal [fr6|-p!q] where semiempirical formula with 
incoherent summation is suggested in application to the DI process. 

In the present paper we address these problems considering HG as the most simple three- 
step event. Being sharply important in itself, it hopefully provides a useful test ground for 
other, more complicated processes. Our study shows that the range of efficiently contributing 
m-channels is rather broad (Sec. |Tp. It is defined by interplay of two factors characterizing 
each m-channel contribution: the well known ATI amplitude (Appendix 0) and much less 
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investigated LAR amplitude (Appendix 0). Importantly, the phases of these amplitudes play 
a major role, showing that the coherence effects are crucial. In particular, solely these effects 
are responsible for the well-known rapid fall-off of HG rates beyond the so called plateau 
domain. At once this finding provokes a new question: how one can simplify summation over 
broad range of essential intermediate m-channels. We achieve this objective (Sec. [TTTD by 
introducing a concept of effective ATI channels, or more briefly, effective channels (EC). Each 
EC is characterized by a complex-valued channel label rric-, i.e. a complex-valued number of 
absorbed photons. The latter is defined from the equation with simple and appealing physical 
meaning. The summation over a large number of intermediate m-channels is replaced by 
summation over very small number of EC (usually one or two). In the higher part of 
harmonic spectrum, beyond cutoff, it is sufficient to take into account a single EC. Two 
EC provide description of the plateau domain, including intricate oscillatory interference 
pattern of HG rates. As shown in Sec. EC approximation ensures good quantitative 
agreement with the harmonic intensities obtained within other approaches |]TT],0,IT3 . 



As the system parameters vary, the EC labels move along characteristic trajectories 
in the complex-m plane. This pattern condensedly expresses theoretical background of HG 
process in various regimes. For instance, transition between the regimes of one and two EC, 
i.e., between the cutoff and plateau domains, corresponds to "collision" of two trajectories in 
the complex plane. The oscillatory behaviour of HG rate in the plateau domain is explained 
by interference of two ECs. Thus ECs provide a unified framework for qualitative assessment 
of three-step processes as well as for their quantitative description. Probably the nearest 
analogy for this theoretical tool comprise the well-known Regge poles |]l9l, i.e., the states 
with complex-valued angular momentum, that give an effective description in the quantum 
scattering theory. 



B. Direct implementation of three-step mechanism in theory 

Now it is worthwhile to look in more detail how the three-step mechanism of HG is 



implemented in the theory of Refs. |]I4|,|T5|. The key formula for the amplitude of A^-th 
harmonic generation, 

4 = 2^4^, (1.1a) 

m 

(K„) , (1.1b) 

presents it as a sum where each term is a product of two amplitudes of physical, fully accom- 
plished and observable processes; no "off-shell" entities appear. The first step is described 
by the first factor ^^^^(Km) which is an amplitude of physical ATI process when after ab- 
sorption of m laser photons the active electron acquires a translational momentum p = K^; 
see detailed description of this amplitude within the Keldysh-type approach in Appendix 
^ The other factor, B ^ m^lo{}^■m) , is a combined amplitude of the second and third steps, 
i.e., propagation and laser assisted recombination (PLAR) amplitude. Under an additional 
approximation it can be factorized into propagation (expansion) factor l/RmfMo describing 
the second step and the amplitude of the third step, LAR, CNmii^m)- 

BNmnoiJ^m) = C]\fm(K.m) ■ (1-2) 
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Rm/io is merely an approximate expression for the distance passed by the active electron 
in course of its laser-induced wiggling motion, see Refs. [ 12, 14, 15 1 and formula ( 3.13c| ) be- 



low. The amplitude CNrni^m) of the physical LAR process describes recombination, i.e., 
transition of electron with momentum from the continuum to bound state. Since the 
continuum state is laser-dressed, the recombining electron can emit the A^-th harmonic pho- 
ton, gaining necessary extra energy from the laser field. The formulas for LAR and PLAR 
amplitudes can be found in the Appendix 

The summation in formula ( |1.1|) runs over a number of photons m absorbed on the first 
step. Thus in the laser-dressed continuum the energy conservation constraint selects the 
discrete set of ATI channels, where the electron has a translational momentum K^. These 
channels serve as intermediate states for the three-step HG process. To specify exactly, 
the absolute value of the electron momentum in m-th ATI channel is defined by the energy 
conservation constraint in ATI as 



Km = V2 {miu -Up + Ea) . (1.3) 

where Ea = is electron energy in the initial bound state. Up = -F^/(4ti;^) is the well- 

known ponderomotive potential, uj = 2'k /T is the laser frequency, T is the period, F is the 
electric vector in the linear-polarized laser wave. ATI can play role of the first stage of 
HG process only if the electron momentum has specific direction, namely is directed 
along F. This ensures eventual electron return to the core that makes the final step, LAR, 



possible as discussed in detail in Refs. [12,14,15|. It is worthwhile to indicate once again 



that our approach presumes the single active electron approximation; the atomic units are 
used throughout the paper unless indicated otherwise. The observable HG rates IZn are 
expressed via the amplitudes as 



7^/v = 



27rc3 



"TV 
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(1.4) 



Vt = Nuo is the frequency of emitted harmonic, c is the velocity of light. 

As already mentioned, the general framework for three-step decomposition of complicated 
laser-induced processes is provided by factorization technique of Ref. |T^. The accuracy of 
this theoretical device is governed by multiphoton nature of the processes: the larger is 
number of laser quanta involved, the more accurate results it provides. To avoid confusion, 
it should be emphasized that the factorization is not related to the conventional perturbation 
theory, or to its simplified version known as the pole approximation. 



Here we only briefly outline some features of the derivation p2| , p!4| , p!5| necessary for 
understanding of the present development. The basic expression for the HG amplitudes 
contains integration over two time variables, t and t' . The factorization technique allowed 
us to disentangle these integrations basing on the adiabatic approach. Mathematically the 
latter implies the saddle point approximation for calculating the integrals which emerge 
when a particular representation for the Green function in terms of intermediate states is 
chosen. The saddle point integration over time variable t' is intrinsic for the factorization 
technique, whereas another time integration could be carried without approximations, for 
instance, numerically. The latter approach was adopted in Refs. P^JT3| and in Sec. |I| of 



the present paper. However, EC approximation implies saddle point integration also over 
time variable t. The details could be found in Appendices ^ and 0. Here we only indicate 
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physical meaning of the emerging saddle points. For t' variable the saddle point t'^^^ , being 



given by Eq. (A7), is the instant of time when the electron is emitted into an intermediate 
ATI channel. The saddle point in t-integral, ti^^^, see Eq. ( [B13| ), is the time of laser-assisted 



recombination of electron with the core, accompanied by emission of high-energy photon. 
Note that generally both t^^^^ and t]s^m^l are complex- valued. 

Formula (|1.1|) is important conceptually being the most direct and fully quantum de- 
scription of the three-step mechanism of HG. It can serve also as a practical computational 
tool. The latter statement was testified by a very good quantitative agreement between cal- 
culations of HG rates by formulae ( |1 . 1| ) , ( |1.4| ) |n,r^| and the benchmark results by Becker 
et al [O] for HG by H~ ion. However, as discussed in Sec. [I A|, it was not explored before 



what is the range of channel label m that gives substantial contribution to the rate. This 
problem is addressed in Sec. |I| where we show that a rather large number of intermedi- 
ate ATI channels is to be taken into account. The effective summation method of Sec. |II 
transforms sum ( |l.la ) over large number of m-terms to the sum over very small number 
(one or two) of EC terms. The illustrative applications of EC approach in Sec. |^ serve 
to demonstrate its reliability. EC representation is convenient also for analysis of general 
features of HG spectrum. The latter is known to consist of three parts. The initial rapid 
decrease is followed by the plateau domain and the rapid cutoff region. Being interested in 
the generation of rather high harmonics, we do not consider below the initial low-A^ part 
of the spectrum. The upper boundary A'^b of the plateau is given by the known expression 



N^uj = -K^ + 3.17f/p . 



:i.5) 



In Sec. |III C| we show how this important result follows from the three-step mechanism cast 
in the EC form. Section ^ contains concluding discussion. 



II. COMPOSITION OF HG AMPLITUDE FROM ATI CONTRIBUTIONS 

Representation ( |1.1|) of HG amplitudes ci^ as a sum over ATI contributions is character- 
ized by a number of interesting features. We illustrate them by Fig. |I| where different terms 
in the sum ( |l.la| ) are quantitatively presented for HG by H^ ion in the laser field with the 



frequency uo = 0.0043 and the intensity / = 10^^ W/cm^. We show separately the squared 
moduli and phases for the ATI amplitude A^^^ (open squares), PLAR amplitude BNm^io 
(open triangles) and the resulting term in the sum ( p.. la]) ci^^ (closed circles) in dependence 
on m for some representative values of the harmonic order N. Namely, A^ = 15, 25 lie in the 
plateau domain of HG spectrum whereas A^ = 39, 51 are in the cutoff region. Only open 



intermediate ATI channels are included in summation ( |l.la| ) that in the present case means 
m > 16. 

The following observations could be made. 



The number of m-terms giving substantial contribution to the sum ( |l.la| ) is essentially 



independent of A^ being about 15 in the example under consideration. 
The domain of m giving substantial contribution is defined, first, by the monotonous 



2 m 

decrease of l^muol with m that ensures effective cut-off from the high-m side. The 
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second factor \BNmfio\ has more comphcated behavior. Indeed, it oscillates with m 
and has substantial magnitude in the much broader range of m than |Am/iol^ (this issue 
is discussed in more detail in the Appendix O). In HG process it is important that for 
high harmonic order N the factor \Bj\fmfio\ is strongly suppressed in the low-m region 
(i.e. near ATI threshold, see Figs. |l]c,d). Therefore, as A^ increases, the domain of 
important contributions to the sum shifts to higher m. 



"'Nm 



decrease with the harmonic 

|2 



Typical values of the modulus squared components 

order A^ increasing. This is explained by decrease in the typical values of the \B]^rn^lo\ 
factor (note that |ylm/iol^ does not depend on A^). Importantly, this relatively slow 
decrease cannot explain the cutoff in the spectrum of modulus squared amplitudes 



N 



For instance, as A^ varies from 39 to 51 the typical values of 



Nm 



decrease 



only by a factor about 3, whereas the resulting rf^ = J2mdtfm decreases by seven 
orders of magnitude ! This clearly indicates that the high-A^ cutoff of HG rate is due 
to the strong cancellation of terms in the sum ( |l.la]) , i.e. that the coherence or phase 
effects are crucial. 

• The plots of the phases show that arg (A^^o) decreases with m whereas arg (-BATm^to) 
increases. The resulting phase arg (ci^vm) = ^'^z{^m^lo) + arg (-BAfm^to) grows with m 
rather rapidly. For larger A^ variation of phases becomes faster. For instance, in the 
case A^ = 51 the phase arg {d^j^ increases by about lOvr along the substantial m- 
domain (20 < m < 35). Namely the rapid phase variation governs cancellations in the 
sum ( p.. lap and rapid decrease of the HG rates in the fall-off region. 

The rapid variation of phases is generally characteristic for the semiclassical dynamics. 
Therefore one can suggest that in our problem some semiclassical-type method of summation 
over ATI contributions could be developed that reflects the physics of the process. This 
program is implemented in the next section. 



III. EFFECTIVE ATI CHANNELS 



A. Poisson summation of ATI channel contributions 



Further in the basic formula ( |1 . 1|) we employ for both ATI and PLAR amplitudes the 
saddle-point approximation given by formulas (|A2|) and (p9|) , respectively. This representa- 
tion possesses important advantage of explicitly exposing large phases of amplitudes which 
are of crucial significance in coherent summation as revealed in previous section. The phase 
factors have form of exponents of classical action for the ATI and LAR processes. By using 
these expressions we present the HG amphtude ( p..l|) as 
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2 EES 



exp 



iS {tmp. 1 t ratio 



(3.1) 



Bearing in mind importance of phase factors we write down them explicitly in Eq. 
introducing 
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s{t,t') = s{t)- s{t') + nt 



dr 



Km H sin u;r 1 — Ea 



(3.2) 



whereas all the rest is collected in the pre-exponential factor Qm/i'- 
'4 1 



2lT 



F ( COS ujt'^^^ - cosut^f. 



X 



■ sin Lut 



(3.3) 



S{t,t') has an appealing meaning of action for HG processes expressed as a sum of actions 
for constituent ATI and LAR processes. 

From the pragmatic point of view the attractive feature of the sum representations ( |1 . 1|) 
or ( p.l| ) is obvious: each term in the sum has a clear and simple analytical expression, that 
helps greatly in numerical calculations. However, since the number of essential terms in the 
sum is large, as discussed in Sec. ||, the behavior of the sum may differ drastically from 
behavior of its individual terms. This makes difficult an interpretation of the results which 
follow from ( |l.lj ) or ( |3.1| ).To overcome this disadvantage it is highly desirable to find some 
other representation for the amplitude which would be devoid of extended summations. One 
can obtain guideline for practical implementation of this idea by recalling the fact that sum- 
mation in ( p. . 1|) or ( p.l|) runs over the spectrum of intermediate ATI states. From the general 
physical principles we know that when the spectrum covers a broad range of energies it 
could be advantageous to work within the time-dependent picture. Mathematically this im- 
plies the Fourier transformation which replaces the spectral quantities mu with appropriate 
intervals of time. 

Having these arguments in mind, we carry out the Poisson summation in ( p.l|) that could 
be looked at as desired "Fourier transformation" 



dt, = 2 I dmY, Qmti exp \-iS (t^^,, t' ) - 2m jm 



(3.4) 



Transition from (|3.1|) to (|3.4|) amounts to replacing summation over the integer m by in- 
tegration over related continuous variable. The price is introducing an extra (and infinite) 
summation over an integer j. The transformation is worthy if the latter sum effectively 
contains less terms than the original m-sum (|l.la| ), i.e., converges more rapidly. 



B. Saddle point method and effective channel representation 

Our next step is evaluation of integrals over m in ( |3.4| ). As demonstrated in Sec. 0, 
the phase of the term c?^^ in the sum ( |l.la|) varies rapidly with the summation index m. 
Physically this follows from a multiphoton nature of the process under consideration. We 
presume that the exponent in (|3.4|) is responsible for this, whereas the pre-exponential factor 
Qmfj, varies slowly. This allows us to use the saddle point method to carry out explicitly the 
integration over the ATI channel label m. The position of the saddle point (s) is governed 
by equation 
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(3.5) 



When taking derivative in Eq. ( p.5[ ) one has to remember that there are two sources of 
m-dependence: first, via t'^^^^ and that are the integration hmits in the definition ( p.2|) , 
and, second, via Km that enters the integrand in 
simphfied by the fact that 



. However the situation is drastically 



(3.6a) 



f'=f' 



d_ 
dt 



Sit,t') 



(3.6b) 



t — tfi 



since t^^^^ and tmij, are the saddle points in integration respectively over t' and t variables as 
discussed briefly in Sec. |I| and in more detail in Appendices ^ and ^ see formulae ( |X5| ) and 
( pil| ). Therefore Eq. ( |3.5| ) takes a compact form {dKm/dm = oj/Km) 



Km {t — t' + jT) H / sin ut dr 

UJ Jt' 

or, more explicitly. 

Km (tmf, - + = ~^ (cos Utmf, - COS tut^ 



(3.7a) 



(3.7b) 



It is to be considered together with formulae ( [A7|) and (P13|) defining positions of saddle 
points in integration over t' and t variables respectively. The unknown variable to be defined 
is the ATI label m(.{N,j). As shown below, only complex- valued solutions are possible. Note 
that m enters ( |3^ ), (|A7|) and ( ^131) only via Km ( 11.31) . Therefore solution of the saddle 
point equation (|3.7b ) amounts to finding complex- valued translational momentum Km^ {N, j) 
of electron in the intermediate ATI channel that gives major contribution to the generation 
of A^-th harmonics. 

Within the saddle point approximation formula ( p.4|) is reduced to 



d 



N 



2 X! X! X! QmcfJ. 
j=—oo rUc jJL 



27r 



exp 



(3.8) 



where 



dm? 



s(t 



t' 



{tr. 



UJ 



K„ 



m=mc 

V2Nuj - K 



T^2 V^cM 
m 



+ 3T) + 



df 

dm 



dt 



IK 



dm 



(3.9) 



m=mc 



The necessary derivatives are straightforwardly derived using formulas (|A7|) and (|B13D 
respectively, as 
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mfio 



dt 



dm 



Km F cos ut: 



dm 



Km F cos ujt 



(3.10) 



Even more appealing form of expression ( |3.8| ) is obtained if one reexpresses the con- 
stituent factors in terms of ATI and LAR amplitudes: 



d 



N 



2 ^ ^ Am,^ fioi'Km^ 



1 



R 



271 



exp{—2TTijmc 



(3.11) 
(3.12) 



Here we omit summation over /i implying that it is absorbed into the sum over m^, since each 
saddle point m-c is obtained for some choice of label /i. The form of expansion (propagation) 
factor 



R 



F (cOSUjtm.t, - COSUt'^^^ 



(3.13a) 



is a natural specification of the more general expression obtained earlier P^JI^ and cited in 
Eq. (p3D. In particular, at the saddle point on can use Eq. ( |3.7b| ) to rewrite propagation 
factor as 

(3.13b) 



R, 



K„ 



(tmcfJ. ^'m.cfJ.0 ~^ jF 



Formula ( |3.11| ) provides the most concentrated expression of the result of present work. 
Apart from the smooth factor S(mc, j), it is fully analogous to our starting point expression 
( |1.1D , but with summation over large number of intermediate ATI channels m replaced 
in Eq. (|3.11| ) by summation over small number of effective channels m^. Note that the 
factorization of propagation and LAR amplitudes, as given by formula (1^), appears as 
an additional approximation in the framework of our previous approach ||l^ , p^ where the 
t-integration was carried out numerically. In the EC formulation the factorization is a 
necessary feature, due to the saddle-point method applied to calculation of integral over 
time variable t, see Appendix p. The form ( p.l3a| ) of the expansion factor is more accurate 
than the approximation 



R. 



F cos ut 



(3.13c) 



m/io 



discussed (within unessential sign) earlier ||14| , [15| and implied in formula (|1.2|) of the present 
paper. 

The transparent interpretation of the saddle point equation ( |3.7a]) is based on the fact 
that 



F /■* 

z (t, t'; p) = pit — t' + jT) H / sin cur dr 

UJ Jt' 



(3.14) 



is the electron displacement along z-axis in course of its wiggling motion in the laser field 
as time varies from t' to t. The axis is directed along the electric field vector F; the electron 
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translational momentum p has the same direction. In terms of z {t, t'; p) ( |3.14]) the saddle 
point equation ( p. 7a ) can be equivalently rewritten as 



+ (' ;A'J=0. (3,15) 



This equation has a lucid physical meaning. In the semiclassical picture of HG, as imple- 
mented by the present theory, the electron emerges from the under barrier at the instant 
of time t^^jj into the m-th channel of ATI continuum and undergoes backward transition 
into the bound state with emission of the fl photon at time tmn- Both these events occur 
relatively close to the atomic core. Eq. ( p.l5| ) is the condition that after propagation in the 
laser field the electron returns exactly [^| to the emission point. It accounts to the fact that 



the propagation time could be augmented by additional j laser field periods T. The integer 
parameter j conjugate to m in the "Fourier transformation" ( |3.4| ) can be named a recur- 
sion number. As discussed above, in fact equation ( |3.15|) defines an appropriate electron 



momentum K^^. The latter proves to be complex- valued as shown below. 

The return condition ( p.l5 ) is universal in the sense that it does not contain explicitly 



the harmonic order (the latter parameter implicitly defines the return time tm/i)- The 
condition of return is intrinsic in the physical picture of the three-step process. Being 
constituent part of the atomic antenna idea of Ref. 0, it appeared in some form in previous 
studies by other authors. The most close links could be established with Ref. 0. This 
paper employs a different representation that does not appeal to the intermediate ATI 
channels. Nevertheless the saddle point method applied in these calculations also leads to 
the equation for the saddle points with the same physical meaning. Moreover, the same 
number of saddle points plays substantial role (namely, one point at the cutoff region and 
two points at the plateau domain, see below). Unfortunately the saddle point evolution with 
varying harmonic order A^ was not described in detail in the cited paper. Note also that the 
cited paper provides the HG rates only within an adjustable normalization factor. 

As already mentioned, the summation over j in the formulas (U), (gj) or ( gTl reflects 



the fact that the time interval between the emission of electron and emission of the high- 
energy photon could be different; namely, several (j) periods of the laser field could be added 
to it. All these events add up coherently as show Eqs. ( |3.8| ) or ( |3.11| ). The shortest time 



interval corresponds to j = 0. Physically one can anticipate that the related contribution 
prevails because of a spread of electron wave packet. Mathematically the term with j = is 
obtained by plain replacement of the summation over m in Eq. (|1 . 1| ) by integration. Further 
on we concentrate on calculation of this term and consider at first (Sec. [Ill (J| ) the related 
saddle points mc{N) = mc{N,j = 0). 

Let us discuss now a restriction on the summation index j in (|3.8| ). In the saddle point 
method one deforms the integration contour so that it passes via the saddle points; only close 
vicinities of the saddle point effectively govern the integral magnitude. Generally, the shifted 
contour passes only via some subset of the saddle points available; only the points from this 
subset contribute to the integral. We deal with multi-dimensional problem (integration over 
variables m, t, t') that casts analysis of the integration contour deformation as a complicated 
task. To circumvent this obstacle we adopt an approach based on simple physical arguments. 
As discussed above, the recursion number j counts a number of periods of time that elapse 
between the electron escape and its return to the atom. The three-step mechanism presumes 
that the moment of escape must precede the moment of return. This condition can only be 
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satisfied for non- negative values of j. Basing on this physical argument we assume that the 
saddle points contribute to the integral only for non-negative j > 0. 

Saddle points which contribute to summation over in formula ( p.8| ) should satisfy 
additional condition 

Im [S (t^^f,, C^^J + 27rjmJ < , (3.16) 

which ensures that the factor in the exponent in the right-hand side of Eq. ( |3.8| ) always 
reduces the absolute value of the amplitude. Again we do not attempt to vahdate ( p.l6| ) 
mathematically, but note instead that this condition is very similar to a conventional re- 
striction on the resonances in stationary processes. Remember that the energy Eq of a 
conventional resonance has the real and imaginary parts. The later one describes the width 
of the resonance ImEo = — r/2 < 0, and must be negative. In the time- dependent formalism 
the resonance contribution to the amplitude is given by the factor exp{—iEQt) , where t > 
is the time elapsed since the resonance intermediate state has been excited. The condition 
Im£^o < selects one of two complex-conjugate poles of the S'-matrix, or Green function, 
thus ensuring the decay of the resonance, i.e., that \exp{—iEot)\ < 0. Similar physical 
meaning has Eq. ( p.l6| ). 



Eqs. ( p.8|) or ( p.ll|) are the major result of this paper. They implement the objective 



formulated in Sec. [Ill A| to carry out the "Fourier transformation" of Eq. ( |3.4| ). The latter 



expression includes summation over the physical, discrete spectrum m of the intermediate 
ATI states. In contrast, Eq. ( |3.8| ) refers to some complex-valued rric which label EC. Addi- 
tionally, Eq. (|3.8|) includes summation over recursion number j, that is a number of laser 
periods T that elapse between the ionization and HG generation. It is important that the 
sum over j converges well. The reason for this originates from the discussed above fact of the 
slow convergence of the spectral representation ( |3.4| ). It is a general, well known fact that 
if the spectral pattern is broad, then the "Fourier transformation" should be well localized. 
Same argument can be presented from another point of view. The shortest (positive) time 
interval between the ionization and HG emission corresponds to j = 0. One can anticipate 
that the related contribution prevails because the spread of the electron wave packet should 
significantly diminish contributions of events with larger j. Presuming that this argument 
is correct, we will replace in the applications (Sec. |IV| ) summation over j in Eq. (|3.11|) by 
the only term with j = 0. Good numerical results obtained illustrate the fact that the term 
j = really dominates. For j = we find only few (one or two) operative EC rric- Thus 
the sum ( |3.11| ) over EC comprises very few terms that is convenient for interpretation of the 
results. 

The physical background of the mathematical transformations above can be briefly sum- 
marized as follows. The electron motion in a laser field in the vicinity of an atom satisfies 
the adiabatic condition because the number of laser quanta absorbed and emitted in the HG 
process is large. The HG amplitude contains a large phase that is identical to the classical 
action S{t,t') (|3.2| ). This phase varies rapidly with all the parameters that govern the elec- 
tron propagation in the intermediate state. Therefore the major contribution to the event 
comes from such situations in which the phase is stationary, i.e. from the saddle points. 
The phase depends on the initial moment of virtual ionization t', the final moment of the 
harmonic emission t and the energy of the electron in the intermediate state. Eqs. ( p.6a| ) 



( p.6b|) , ( |3.5|) represent the saddle-point conditions over these three variables. Altogether 
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they define the two moments of time and the electron energy in the intermediate state. 



C. Analysis of saddle point equation 



By summing and subtracting Eqs. ( [AT] ) and ( [B13D we obtain 



sin utmf, + sintut'^^^ + 2Km ^ = 2Z , 
sin ujtmf^ - sin ut'^^^ = 2Z' 



with 



UJ 

F 



Z' = - (^/2Nuj - - IK 



UJ 

F 



If 2A^a; > k^, then one has Z* = Z'. It is convenient to switch from t'^^^^ and tm^i to 



X 



1 / 



1 / 
-Uj(t 

9. V 



(3.17a) 
(3.17b) 

(3.18) 
(3.19) 



and present ( p.l7a| ) and (|3.17b| ) in the compact form 

UJ 



sin X cos y + K„ 



Z 



cos X siny = Z' . 

In the same notation Eq. ( p.7b| ) reads 

sinx siny 
F y + jrc 

By excluding Km we obtain a set of two equations for two variables x and y: 

' siny ^ 



(3.20a) 
(3.20b) 



(3.21) 



sm X cos y 



= Z , 
y + jTc J 

cos X sin y = Z' 



(3.22a) 
(3.22b) 



The next step is to exclude x that after some algebra gives a compact transcendental equation 
for a single variable y = y + jir: 



cot y — — ] sin^ y — [Z'^ 



(3.23) 



As soon as y is found, one obtains K^, from 



^Km = - 

F 1 — 2/ cot y 



(3.24) 



The variable y (|3.19| ) has a very lucid meaning: 2yjuj is the active electron excursion time 
in the continuum, i.e., the time interval between the first step in the three-step process, ATI 
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and the last step, LAR. Change in the recursion number j is equivalent to change of the 
excursion time {tmn — ^m^o) integer multiple of laser period T. To simplify subsequent 

analysis we restrict it to the case j = 0, or < y = y < ir, that gives the major contribution 
to the amplitude, as discussed in Sec. [1IIB| . 



When generation of high harmonics is considered, 2Nlo ^ k^, the estimate 

V2Nu - /«2 > K (3.25) 



is valid that implies that one can put Z' ^ Z ^ V2Nuj - K^uj/{2F). This allows us to 
rewrite equation (|3.23| ) governing saddle point positions in an approximate form 



f(y") = (2^^ - ' (3-26) 
with the universal function /(y^) of a dimensionless parameter y defined as 

(vcotv — 1)^ sin^V 

f{y ) = — — - — — ■ (3.27) 

(ycoty-l) +2/2 

A plot of this function (Fig. |^) shows that in the interval of interest, < ^ = < tt^ it 
has a single maximum located at = 4.173 {i/m = 2.043). For sufficiently small A^, when 
the right hand side of Eq. ( |3.26| ) is less than Ci = f{C,m) = 0.3966, this equation has two 
solutions designated as A and B in Fig. ^ where this situation is explicitly displayed. One 
can verify that these solutions satisfy condition ( |3.16D and therefore both of them contribute 
to the HG amplitude. Note that both solutions correspond to real but different excursion 
times 2y/u!. 

When A^ increases, the roots A and B come closer and eventually merge at some critical 
value of A^ = A'b defined from 



2 



4F2 



2Ni,LU - K^) = Ci , (3.28) 



that is 



1 2F'^ 1 

Nkuj = -n^ + Ci^ = -K^ + 3.1731C/p . (3.29) 

2 iu 2 

For larger N > Ny, the solutions A and B spit again, but this time they result in complex- 
valued y. The way to see this, is to note that in the vicinity of the maximum the function 
/(?/2) behaves as /(y^) — fiVm) ~ ~iy ~ VmY- Further straightforward analyses indicates 
that only one of these two complex- valued solutions satisfies condition (|3.16|) and contributes 
to the amplitude, while another one gives no contribution. 

Thus the number and the character of relevant solutions of Eq. ( p. 261) differs for N < 



and A^ > A^b- The excursion time switches from real to complex-valued that indicates that 
the HG process changes it nature from the classically allowed to classically forbidden one. 
Accordingly, as detailed in Sec. 0, at the point N < N\y the HG spectrum undergoes a 
transition from the plateau to the cutoff region. Thus Eq. ( |3.29 ) rederives the well-known 



upper border of the plateau domain in the HG spectrum ( |1.5| ) with a slightly different 
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coefficient in front of Up. The validity condition for this derivation (|3.25|) essentially means 
that the photon energy A^b"^ is niuch higher than the initial electron binding energy 

Before concluding this subsection we remark that one more solution of Eq. ( p.26| ) is shown 
in Fig. ^ being labeled as C. It corresponds to negative y"^ and has simple A^-dependence 
since the function /(y^) is monotonous on the semiaxis < 0. The fact that the solution 
C results in imaginary excursion time 2y/uj indicates that it gives small contribution to the 
HG amplitude. 



IV. ILLUSTRATIVE APPLICATIONS 

A. Scanning harmonic order for fixed laser intensity 

There are two natural and complementary outlooks on the results for HG process. One 
can fix the laser intensity / and scan the rates for harmonics of different order, as done in 
this subsection, or one can fix the harmonic order N and vary the laser intensity (Sec. [IV B| ). 



In both cases we consider HG by H~ ion in the laser field with the frequency u = 0.0043. 

Fig. D depicts a typical pattern of evolution of mdN) in the complex-m plane for fixed 
laser field intensity / = 10^^ W/cm^. Positions of three complex-valued roots of the exact 
saddle-point equation (|3.7b|) are shown by symbols of different shape. The roots move as 
the harmonic order varies; we show their positions for odd (physical) values of A^. For 
large A^ only one root has negative imaginary part (it is shown by diamonds in Fig. ^ 
that ensures harmonic rate decrease with increasing A^. Fig. ^ shows results of the rates 
calculation in the saddle point approximation ( p.8|) where only this single root is taken 
into account in the summation over rric- The cutoff region in the rates spectrum is nicely 
reproduced, as well as the overall pattern in the plateau domain. However the single- saddle- 
point approximation does not reproduce some structures in the A^-dependence of HG rates. 
In the plateau domain the saddle point m(.{N) moves close to the real m-axis. For small A^ 
this saddle point approaches the point m = rrith on the real axis in the complex-m plane 
that corresponds to the ATI threshold 

mth = — + ^ . 4.1 

Zuj 00 

Transition to the cutoff region with A^ increasing corresponds to steep bending of the 
saddle point trajectory after which it moves almost perpendicular to the real axis. The 
reason of the trajectory bending is seen to be a "collision" with another saddle point shown 
by the circles in Fig. |[ In the "collision" region the two saddle points ("diamonds" and 
"circles" ) abruptly change direction of motion forming the characteristic cross-like pattern. 
The "collision" occurs at A^ = A'coi ~ 35. Approximately one can identify A"coi with A^b- For 
A^ < A^coi the two saddle points rric are to be taken into account in formula (|3.8| ). As shown in 
Fig. ^, this improves the results for the rates in the plateau domain by producing appropriate 
structures in the A^-dependence. Remarkably, such a two-saddle-point calculation gives 
correct positions of minima and maxima in the rate A^-dependence, albeit the magnitudes of 
the rate variation is reproduced somewhat worse; for instance the depth of the minimum at 
A^ = 17 is quite strongly overestimated. Tentatively one can attribute this to the fact that 
if the principal term in the approximation for considered here proves to be anomalously 
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small for some A^, then the correction terms omitted in our calculations become relatively 
important. The good overall description of the structure unambiguously identifies its nature 
as a result of an interference between the contributions coming from two effective ATI 
channels. 

The plain saddle point approximation ( p.8| ) presumes that the saddle points are well 
separated from each other; otherwise the more complicated uniform approximations are to 
be constructed. In our problem this refers to the "collision" region ^ A"coi where our 
simple approximation somewhat overestimate rates; however we do not resort here to the 
more sophisticated mathematical constructions. Note that the third saddle point shown by 
triangles in Fig. ^ is substantial only in calculations for the lowest harmonics. Therefore it 
is not taken into account in the present calculations. 



B. Scanning laser intensity for fixed harmonic order 

It is equally instructive to see how the rate of some individual harmonic depends on the 
laser intensity /. Fig. ^ shows trajectories of the saddle points as / varies; only two most 
important saddle points are considered. For small intensities only the saddle point shown 
by diamonds is operative in evaluation of HG rate. This corresponds to the rapid-fall-off 
regime beyond the plateau domain. Here Im is large and negative. As / increases, the 
saddle point moves to the real m axis. After crossing the axis the trajectory abruptly bends 
at / about 3 x 10^° W/cm^. For higher intensities, in the plateau domain, the second saddle 
point, shown by circles, should also be taken into account in calculation of HG rates. In this 
region the single-saddle-point approximation correctly reproduces an average behavior of the 
HG rate, whereas the two-saddle-point approximation shows also characteristic oscillations 
of the rate with / (Fig.|^). The positions of minima and maxima are correctly reproduced 
for medium intensities. On the higher I side an additional maximum at / = 9 x lO^'' W/cm^ 
is beyond the two-saddle-point approximation being tentatively due to the contribution of 
the third saddle point omitted in the present calculations. This contribution could be also 
the source of an additional minor structure in the medium intensity region. It could explain 
also the fact that the two-saddle-point approximation gives deeper minima in the rates than 
the numerical summation in ( [1 . 1|) . 

The present calculations clearly show that the structures in HG rates as function of / 
is related to the interference effects, in agreement with conclusion reached by Lewenstein 
et a/ 0. The alternative explanation |TO|,|TT| relates this structure to the threshold effects, 



namely to the successive closure of ATI channels by increasing ponderomotive potential as 
I grows. However, the threshold effects should be manifested also in the single-saddle-point 
approximation; Fig. ^ demonstrates that this is not the case. 

In the plateau regime the saddle point rric shown by diamonds in Fig. ^ lies in the 
complex-m plane close to the threshold value mth calculated for the same laser intensity 
using formula ([4.1| ). As I increases, tends to approach mth- This behavior is in agreement 
with the three-step model developed by Corkum who presumed that the ATI electron 
leaves an atom with zero velocity. Note however, that this presumption is valid only within 
single-saddle-point approximation (shown by diamonds in Fig. |^) and only in the plateau 
domain. The point m = mth is a branch point for any function which depends on m 
via Km- In particular, generally the preexponential factor Qmfi (^) in expression ([ 
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has such a branch point. In principle one can look for some refined version of the saddle 
point approximation which accounts for closeness of a saddle point and a branch point. 
Note however that such modifications infiuence only the preexponential factor but not the 
exponent which is the principal object of interest. As discussed above, description of the 
transition between one-saddle point and two-saddle point regimes is more difficult for the 
semiclassical-type theory since it requires more sophisticated approach. In the present case 
the situation is aggravated by the presence of the branch point. It makes behavior of 
trajectories near the 'collision point' quite different from the conventional cross-like pattern 
(Fig. ^) and leads to appearance of a spike on the dashed curve in Fig. ^ at the borderline 
intensity / = 2.7 x 10^° W/cm^. 



V. CONCLUSION 

The present study has an objective to get a better insight into the three-step mechanism 
of processes in strong laser field. We take HG as the simplest case of three-step process and 
employ representation of its amplitude in terms of amplitudes of physical, fully accomplished 
ATI and PLAR processes. The number of alternative paths in the three-step picture is 
generally infinite, each path being labeled by the number m of laser photons absorbed in the 
first step, ATI. We explore the range of substantial intermediate ATI channels and reveal 
crucial role of coherent interference. The effective scheme of summation of the contributions 
coming from different intermediate ATI channels is developed based on the saddle point 
method. Physically this approach is justified by the multiphoton nature of the process. 
Due to it in the strong-field regime both ATI and PLAR amplitudes have phases that vary 
rapidly with m. Within this framework we develop concept of effective ATI channels that 
is important basically and useful for practical calculations. In this approach interference of 
infinite number of competing paths is replaced by the single-path picture for generation of 
high harmonic, or by interference of two paths for harmonics lying in the plateau domain. 
In particular, the structure in HG rates in the latter domain is understood as the simple 
two-paths interference pattern. 

The effective ATI channels are related to the particular classical electron trajectories that 
ensure that electron emitted from the atom at the first, ATI stage of HG process returns to 
the core to make possible subsequent laser assisted recombination. The possible trajectories 
possessing this property in principle comprise an infinite set differing by the excursion time 
in continuum which is integer multiple of laser period T. However due to quantum wave 
packet spread only the trajectories with the shortest excursion time provide substantial 
contribution to HG rate. The crucial step in this development is compexification of the 
problem: the initial and final time is complex-valued that makes entire classical trajectory 
also complex-valued. The simplicity of description achieved along these lines is in sharp 
contrast with the approach based on the conventional real-time classical trajectories where 
the problems related to chaotic irregular classical motion emerge in full scale ^3 . 



16 



ACKNOWLEDGMENTS 



This work has been supported by the Austrahan Research CounciL V. N. O. acknowl- 
edges the hospitahty of the staff of the School of Physics of UNSW where this work has been 
carried out. 



APPENDIX A: ATI AMPLITUDE IN ADIABATIC APPROXIMATION 



The original expression for the ATI amplitude within the Keldysh |jT]] approximation 
reads 



T 



dt' ($p(t') I VF{t)de I ^ait')) 



(Al) 



where $p is the Volkov state of the electron with the translational momentum p, $a is the 
electron initial bound state. Details of all the definitions could be found in Ref. iTHl; here we 



indicate only that Vp(r,t) = r ■ F cos cut describes an interaction of the active electron and 
the laser field, with F and T being introduced in Sec. [TB] . The absolute value of the 
electron momentum in the m-the ATI channel is subject to energy conservation constraint 
being given by formula ( |1.3|) {pm = Km)- 



The adiabatic treatment of the ATI process was developed in Ref. and subsequently 
applied in Refs. p2| , p!4| , |l5i . It presents the ATI amplitude A^(p) as 



^K(p) 



[270 

T 



exp 



(A2) 
(A3) 



where S{t) is the classical action 

1 /•* / F \^ 
5'(t) = -y (ir (p + — sin cur 1 -Eat. 



The position of the saddle points t'^^ in the complex f plane is defined by equation 
or, more explicitly. 



p+-sinw4^) +K 



UJ 







(A4) 



(A5) 



(A6) 



Expression for HG amplitude ( |1.1| ) refers to the particular saddle point (/i = /io) defined for 
the monochromatic laser field by the expressions 
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sin MO = -J ("^™ + ' (^''^) 



2 



cosw4^g = + . (A7b) 

The other notations are as follows: is initial electron binding energy, i/ = Z / Z \s 
the charge of the atomic residual core (z/ = Z = for a negative ion), /, are the active 
electron orbital momentum and its projection in the initial state. In the present context 
the unit vector p coincides with F = Y/F. The coefficients Aa are specified in Ref. ||21| 
being tabulated for many atoms and ions |2^. Mathematically formula (|A2|) is obtained 
by using the saddle point method to carry out integration over time variable t' in formula 
The saddle point positions are defined by Eqs. ( |A5|) or ( |X6| ) that are to be considered 
together with the energy conservation constraint ( |1.3| ). The summation in ( [A^ ) runs over 
the saddle points t'^^ in the plane of the complex- valued time t' . The saddle points t',^^ lie 
symmetrically with respect to the real t' axis. For the monochromatic laser field there are 
four saddle points in the interval < Ret^^ < T, two of them lying in the upper half plane 
(Imt^^ > 0). Only these two saddle points are included into the summation in (|A2|) . If ATI 



is considered as the first step in HG process, then, as discussed in detail in Ref. ||T5|, only 
one of these two saddle points is effectively operative, namely, that specified by formulae 



APPENDIX B: LAR AND PLAR AMPLITUDES IN ADIABATIC 

APPROXIMATION 



The PLAR amplitude defined as 



dt 



{^a{t) I exp(^^]^)4 I $p(t)) 



differs from the LAR (Laser Assisted Recombination) amplitude Cm(p) 



(Bl) 



2-kT 



dt{%{t) I exp(2f]t)4 I <^a{t)) 



(B2) 



only by the "propagation" or "expansion" factor 1/ Ro{t,t'^^g) in the integrand. The 
distance passed by electron in the laser field between initial time t' and final time t is 



approximated as [15 



Roit,t') 



(cos Ojt — COS Ojt') 



(B3) 



In (pip, ( [B2[ ) (ig = e ■ r is the dipole interaction operator. 

Consider electron in the laser-dressed continuum state with the translational momentum 
p. It can recombine to the bound state with emission of the photon. Possible frequencies 
of the emitted photon form an equidistant pattern: 
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^^m'' = ^p' + f^p-^a + Mo; (B4a) 



with an integer M. Expression ( |B4a| ) can be reparametrized to the form 



n)^^={J + T^)uj , (B4b) 

where J is another integer and the fractional parameter 77 (0 < 77 < 1) is governed by the 
value of the initial momentum p. In the zero-laser-field limit [F — > 0) only emission of the 
photon with the frequency f^p'^ = | — i^^ is allowed. The presence of intensive laser field 
makes possible the processes when laser photons are absorbed from the field or transmitted 



to it, and thus the entire spectrum (|B4|) is produced. Such a LAR process has not yet 
received much attention in literature; as far as we know, our recent study p5| is the only 
theoretical paper on the subject. In it the initial electron translational momentum p was 
considered as an arbitrary input parameter as required in applications to LAR process in 
laser plasma. In other terms, the fractional parameter r] was arbitrary. 

Here we are interested in LAR as a constituent part of HG process. This application 
has some particular features. First, all the frequencies in the photon spectrum are integer 
multiples of uo, i.e. 77 = 0. This happens because the initial translational momentum p of 
LAR process could not be arbitrary since only the discrete subset of continuum states is 
populated by ATI from the bound state. Namely, the translational momentum is subject to 
the constraint ( |1.3| ) with m being the number of laser photons absorbed on the initial ATI 
stage of HG process. 

Second, when an individual harmonic is considered, its order is fixed, and the pa- 
rameter m is scanned when summation is carried out according to the expression ( |l.laD . 



Equivalently, the electron momentum p = Km, is scanned along the discrete set of allowed 
values Km, see Eq. (|1.3| ). This is in variance with the laser plasma applications when it is 
natural to presume that p is fixed and Q^^^ is scanned. The third special feature is that 
the electron translational momentum p = is parallel to the electric field amplitude F. 

Finally, the object of interest in HG theory is the PLAR amplitude ( [BT| ) rather than 
LAR amplitude (p2D, although both amplitudes are quite similar. Below we briefly expose 
modification of some results in the LAR theory relevant to the present application. 

The expressions ( [BTD and (|B^), respectively, for PLAR and LAR process amplitudes 
are valid in the Keldysh-type approximation. By using the Fourier transformation they are 
rewritten as 

T 

B^m^oi^m) = -7^ I dt \^ exp{7[^]t-g(t)]} (-K^--sina;t) , (B5) 

CNmiKm) = - ^ Jdtexp {i [nt - S{t)]} (-Km - ^ sinujt] , (B6) 
^ ^ 

where the classical action S{t) is introduced above by formula (|A4D . The function (j)^\q) is 
defined as 

0i^)(q)=7(e-Vj0.(q) . (B7) 
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where 0a (q) is the Fourier transform of the bound state wave function </)a(r): 

0a(q) = J exp(-iqr) ^^(r) . (B8) 

The time integrals in ( [B5| ) and ( PS| ) could be evaluated using the saddle point approxi- 
mation: 



1 Lu"^ ~ / F ^ 

BNm^O = ^ E —f ~^a^ " - siu O^t^^ 

^ttJ F (cost^t;„^j, - coscjt^^j V ^ y 



X 



X y exp {z [iVa;t„^ - S{t^^)]} . (B9) 



Cjv^i = E <^a^ ^-Km - ^sintut^^ 



X 



X ) exp {z [iV^t„^ - S{t^^)]} , (BIO) 

where summation is to be taken over the saddle points tmii operative in the contour integra- 
tion. The position of saddle points in the complex t-plane is governed by the equation 



or, more explicitly 



^'(t„^)-fi = , (Bll) 



^ fp+ ^sincjt^^ ) =Ea + n. (B12) 



sin cjt„^ = ^ (-Krn ± ^2Nuj - H?) . (B13) 




Its solution is 



For = the formula ( P13|) coincides with that governing the saddle point position in 
case of ATI processes [see Eq. ([A7|)]. As in the latter case, Eq. (pil|) has four solution per 
the field cycle (i.e for < Ret^^ < T). Physically t^^^^ is the time when the electron is 
emitted at the first step of HG process and tmij. is the time of electron return back to the 
initial bound state at the last stage. Note however that both t'^^^ and tm^t are generally 
complex- valued. 

Formula ( |B13| ) might give real- valued saddle points that correspond to the classically 



allowed LAR or PLAR. According to Eq. (|B12|) this means emission of the high-energy 



photon at the real moments of time when the instantaneous kinetic energy of the electron 
in the laser field differs from the bound state energy Ea by Q. In this respect LAR or PLAR 
processes differ basically from ATI process which is always described by the complex-valued 
saddle points that corresponds to the classically forbidden or tunneling transitions. The 
classically allowed transitions in LAR or PLAR are absent for small N when < /t^/(2c<j). 
However we are interested only in generation of harmonics with sufficiently high order: 



N > — . (B14) 
2uj ^ ' 
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In this case the classically allowed population is operative provided m is not too large, 
namely 



Km<- + V2Nuj - k2 , (B15a) 

CO 



m< m% = N + 3^ + — V2Nuj - . (B15b) 



I.e. 

If is sufficiently large, namely if 

F 



— < V2Nlu - k2 , (B16a) 

UJ 



I.e. 

k2 2a 



N>— + ^, (B16b) 

ZuJ UJ 

then the domain of classically populated m-channels is bounded from below by the condition 



Km> + ^2Nuj - , (B17a) 

UJ 



I.e. 



m>m^ = iV+^- — V2mJ^^ . (B17b) 

UJ UJ^ 

In the interval < m < the classical population is governed by two real saddle points 
tm/j. per field cycle. Note that for the real LAR or PLAR process one always implies that 
the m-th ATI channel is open, i.e. m is sufficiently large: 

m > mth (B18) 

where mth is defined by Eq. ( [4.11 ). The same constraint is assumed also in the sum over m 
in formula (|l.la|) implementing the three-step mechanism of HG 0,115 . 
For intermediate N lying in the interval 

!^<N<— + ^ (B19) 
2uj 2uj UJ 

another regime of classical population becomes possible. Namely, if the condition 

F 



I.e. 



K„,< V2Nuj - /t2 , (B20a) 

UJ 



OTJ p 

mtY,<m<m^ = N+—^ V2Nuj - (B20b) 

UJ UJ^ 



is satisfied, there are four saddle points t^u per field cycle. 
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In Fig. ^ we illustrate these results taking LAR of electron on hydrogen atom with 
formation of negative H~ ion. The laser field intensity is J = 10^^ W/cm^. Open ATI 
channels lie at m > 16. The intermediate regime ( |B19| ) takes place for 7 < N < 24. In Fig. 
^ we show m distributions of LAR modulus squared amplitudes |C7VmP for two adjacent 
A^ from the intermediate- A^ domain, A^ = 13 and A^ = 15. For A^ = 13 the classically 
forbidden domain is m > 62; the classical population in the two-saddle-point regime is 
operative for 19 < m < 61, the four-saddle-point regime is operative for 16 < m < 18. For 
A^ = 15 the classical population is forbidden for m > 66; it is allowed in the two-saddle-point 
regime for 18 < m < 66; the four-saddle-point regime is operative for 16 < m < 17. The 
region of transition from classically allowed to classically forbidden LAR requires some more 
elaborate treatment to reproduce related transition pattern similar to Airy function. The 
special analysis is required also to describe a transition from two-saddle-point to four-saddle 
point regime. However here we do not pursue the objective of detailed uniform description 
of LAR or PLAR amplitudes. 

Note that for A^ = the saddle points coincide with the singularity of the integrand 
(l)^\q). This coincidence occurs also in the adiabatic treatment of ATI. The difference is 
that in the particular case of above threshold detachment of negative ion the singularity 
in the integrand is canceled and the standard form of the saddle point approximation is 
applicable. Otherwise (i.e. in case of non-zero atomic core charge) the saddle point method 
is to be applied in the somewhat modified form. In case of LAR the singularity is always 
present, and hence the modified form of the saddle point method should be applied for 
A^ = 0. Respectively, the small A^ case requires some special treatment (some sort of 
uniform approximation; note that all these modifications affect only the pre-exponential 
factor, but the principal exponent remains the same). Below we are interested in the large 
A^ case that allows us to put aside these complications. However, it is to be remembered that 
coincidence of the saddle point and the singularity of the momentum space wave function 
means that only large-r asymptote of the initial state wave function V5a(r) in the coordinate 
space is of importance. This is the case for ATI, but not for LAR or PLAR, as follows from 
the preceding discussion. In other words, the LAR amplitude (or its PLAR counterpart) 
generally are more sensitive for the short-range behavior of the wave function. 
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FIGURES 



FIG. 1. Contributions to the HG amplitude d'^ from different ATI m-channels according to 
the expression (1.1). The H~ ion in the laser field with the frequency w = 0.0043 and intensity 
/ = 10^^ W/cm^ produces harmonic of the order N indicated in the plots. The square moduli 
|MmP and phases arg(Mm) of three different amplitudes are pictured. The closed circles 
show the squared moduli and phases of the terms = Amiif) -BiVm/io in the sum (LI) over ATI 
channel label m. We present also the squared moduli and phases of the physically meaningful 
factors that constitute d~^^: ATI amplitude Am^g (open squares) and PLAR amplitude BNmtJ.Q 
(open triangles). For graphical representation the convenient scaling factors are introduced in the 



plots of squared moduli, namely as |MmP we show everywhere 10^ 



- 10^ l^mMoP and 10^ \B 
10"^ l^m^ol^ and 10^ \B 



N miiQ I j 



for N = 25 



N m^Q I 



for iV = 51 - 3 



5 • 10^ 

10^ 



2d+ 



and also: for = 15 



Am ^lol'^ and 10"^ I^Arm^iol^; for A^ = 39 
^m/xol and 10^ \BNm^io\ ■ 



FIG. 2. The universal function /(^ (|3^ . 



FIG. 3. Trajectories of the saddle points mc{N) in the complex-m plane for fixed laser 
intensity / = 10^^ W/cm^ and varying harmonic order A^. The results are for H" ion in the laser 
field with the frequency lv = 0.0043. Positions of three saddle points for odd integer A^ are denoted 
respectively by diamonds, circles and triangles. The plot (a) gives general overview, the plot (b) 
shows the region where two saddle points "collide", whereas the plot (c) details behavior of the 
saddle point trajectory in the vicinity of m = 772-th- 



FIG. 4. Harmonic generation rates ( |1.4| ) (in sec ^) for H ion in the laser field with the 
frequency u = 0.0043 and various values of intensity / as indicated in the plots. Closed circles - 



results obtained by Becker et al |Tl|, open circles - our calculations |15| in the dipole-length gauge 
performing numerical summation ( |1 . ![ ) over contributions of different ATI channels, open diamonds 
- present results within the saddle point approximation (|3.8| ) for the summation over ATI channels 
contributions with a single saddle point mc{N) taken into account (namely, the saddle point shown 
by diamonds in Fig. ^); open squares - same but taking into account two saddle points mc{N) 
(namely, these shown by diamonds and circles in Fig. ^). 



FIG. 5. Trajectories of the saddle points mc{N) in the complex-m plane for fixed harmonic 
order A^ = 15 and varying laser intensity I = ( x 10^^ W/cm^ (HG by H~ ion in the laser field 
with the frequency u) = 0.0043 is considered as in Fig. |3|). The plot (a) gives general overview, 
the plot (b) shows vicinity of the real-m axis. Positions of two saddle points for odd integer N 
are denoted respectively by diamonds and circles with the numbers indicating value of the factor 
i.e. intensity in the units 10^*^ W/cm^. The crosses in the plot (b) show positions of the ATI 
threshold defined by formula ( [4.1[ ) for the same values of 
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FIG. 6. The rate parameter M^r = 2 log 



10 



4 



for the harmonic = 15 as a function of 



the laser field intensity I (HG by H ion in the laser field with the frequency lv = 0.0043). Solid 

over 



curve - our calculations in the dipole- length gauge performing numerical summation (1.1 



contributions of different ATI channels, dashed curve ~ present results within the saddle point 
approximation (^^) for the summation over ATI channels contributions with a single saddle point 
mc{N) taken into account (namely, the saddle point shown by diamonds in Fig. |5|); dotted curve 
- same but taking into account two saddle points mc{N) (namely, these shown by diamonds and 
circles in Fig. |3|). The bars with numbers m indicate the threshold intensities Im such that for 
I > Im the m-th ATI channel is closed due to ponderomotive potential. 



FIG. 7. Squared modulus of the amplitude Cjvm of the laser assisted recombination of electron 
into bound state in the laser field with the frequency uj = 0.0043 and intensity / = 10^^ W/cm^. 
The harmonic order is fixed [N = 13 (circles) and = 15 (triangles) in (a) and = 39 (circles) 
in (b)]. Closed symbols - calculations using the formula ([B6|) with numerical integration over t; 



open symbols - adiabatic approximation ( BIO ) with various number of saddle points taken into 
account for different m as described in the text. 
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